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Abstract. 

We develop an analytical theory for diffusive random lasers by coupling 
the transport theory of the disordered medium to the semiclassical laser rate 
equations, accounting for (coherent) stimulated and (incoherent) spontaneous 
emission. From the causality of wave propagation in an amplifying, diffusive 
medium we derive a novel length scale which we identify with the average mode 
radius of the lasing quasi-modes. We show further that loss at the surface of 
the laser-active medium is crucial in order to stabilize a stationary lasing state. 
The solution of the transport theory of random lasers for a layer geometry 
with appropriate surface boundary conditions yields the spatial profile of the 
light intensity and of the population inversion. The dependence of the intensity 
correlation length on the pump rate is in qualitative agreement with experimental 
and numerical findings. 



1. Introduction 

A random laser is a system formed by randomly distributed scatterers embedded 
in a host medium where either scatterer or host medium or both provide optical 
gain through stimulated emission [T] . Recent observations of random lasing in a wide 
variety of systems, like powders of semiconductor nanoparticles [2[[3l[4l[5], organic dyes 
in strongly scattering media [6j E3 IS] , organic films or nanofibers (SUTUlllI] and ceramics 
[12] . have triggered a rapidly growing interest. For reviews with comprehensive lists 
of references see p~3l [14] . Random lasers share some properties with conventional 
lasers, like threshold behavior [3], narrow spectral lines [TS], or photon statistics 
[TBI [17] , but also exhibit distinctly different properties like multidirectional emission. 
Coherent feedback has unambiguously been demonstrated to be present in strongly 
disordered random lasers [16]. It requires the light to be sufficiently confined in the 
random system. While spatially confined regions from which the laser emission takes 
place have been observed experimentally [TS] [T3], the physical origin of coherent 
feedback, of localized quasimodes and the dependence of their size on the pump 
rate have remained controversial. The possible theoretical explanations range from 
preformed random microresonators |18) to multiple random scattering (diffusion) [19!, 
possibly enhanced by self-interference [15] of waves and the resulting onset of Anderson 
localization (AL) [20]. Conversely, it is an interesting fundamental question how AL 
of light, which has been understood[2TJ [22] as a consequence of self-interference, is 
influenced by the coherent, stimulated amplification in the lasing state. Problem of 
the intensity distribution in a diffusive random laser has been attacked theoretically 
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by phenomenolocical diffusion modeis [23] and numerical calculations in one [2U [25] 
and two spatial dimensions [STJ 12H1 EH]- However, the experimentally observed 
decrease of the lasing spot size with increasing pump rate, has not been explained so 
far. 

In this paper, we address the question of the size of lasing spots in diffusive 
random lasers by an analytical transport theory. We begin the analysis in section 2 by 
an outline of the transport theory of light in a disordered medium with linear gain (i.e., 
constant-in-time amplification rate), including self-interference (so-called Cooperon) 
contributions. By a phenomenological analysis we show that causality implies in the 
presence of linear gain a novel length scale which is to be identified by the average 
radius (spot size) of a lasing mode, R s . In section 3 we extend the theory for a system 
with linear gain to a transport theory including non-linear gain due to stationary 
lasing above threshold. Observing that a stationary lasing state is possible only if 
the amplification in the medium is compensated by loss at the surface of the system, 
we consider a model for a finite-size random laser with infinite extension in the (x,y) 
plane, but finite, constant thickness in the z direction, a geometry relevant for many 
experimental systems [3j [H [HI [16] . Coupling the diffusive transport theory to the rate 
equations of a four-level laser in the stationary state we derive an analytical expression 
for the intensity correlation length £ to be identified with the average lasing spot size, 
R s . Due to the surface boundary conditions the spatial extention of a lasing spot in 
the (x,y) plane obtains a z-dependent profile , £(£). We also analyse its dependence 
on the pump rate. The conclusions are drawn in section 4. on the pump rate and on 
the depth of the mode along the z direction. 



2. Transport Theory for a Diffusive, Linear-gain Medium and Causality 

The propagation of light is described by its wave equation. Neglecting the polarization 
degree of freedom we consider in the following the scalar wave equation for the field 
ty. It reads, 

%r e(f )V U (? ) + V 2 * w (f ) = -tu^j u (f) , (1) 
c c 

where c denotes the vacuum speed of light and ju>{?) an external current source. The 

dielectric constant is e(r)=efc + AeV(r), where the dielectric contrast between the 

background, efc, and the scatterers, e s , has been defined as Ae = e s — e^. The spatial 

arrangement of the scatterers is described through the function V (r)=^2 g (f—R), 

with Sj^ (r ) a localized shape function at random locations R. Linear gain (absorption) 

is described by a temporally constant, negative (positive) imaginary part of and/or 

e s . 

In Refs. [3S[ [301 [31] we have developed a theory for light transport in disordered 
media with linear gain or absorption. It results in an energy-density correlation 
function P£(r — r 1 ^ — t'), which describes how the energy density of the light field 
with frequency u> propagates diffusively between two points in space and time, (r,i), 
(r*, f'). The Fourier transform of the energy- density correlation function P£(Q, O) is 
obtained as 

P£(Q, n) = 3-^ s- ( 2 ) 

where the expression for the coefficient Np is given explicitly in Ref. [3D], but is not 
relevant for the present purpose. The denominator of Eq. ([2]) exhibits the expected 
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diffusion pole structure with the diffusion coefficient D. In addition, in the case 
of a non-conserving medium, i.e. net absorption (gain), there appears the (purely 
imaginary) term i"f a = i^~ 2 D, which has a positive (negative) imaginary part and does 
not vanish in the hydrodynamic limit, f2 — > 0, Q — > 0. The self-consistent solution of 
the transport theory including self-interference of waves (Cooperon contributions) (see 
Refs. [3(31 El]) shows that in the presence of absorption or gain the diffusion coefficient 
D cannot vanish and is in general complex. Hence, truly Anderson localized modes 
do not exist in this case. 

For the case of absorption (j a > 0) it is seen by Fourier transforming Eq. @ 
w.r.t. time, P£(Q,t) = iNpe~(® D +7o)* j that Re"f a represents the loss rate of the 
photonic energy density due to absorption in the medium. Fourier transforming, on 
the other hand, Eq. @ w.r.t. space in the stationary limit (SI —> 0)) shows that 
£ a = Re\J (jjD) is the length scale over which the energy density of diffusive modes 
is correlated in the lossy medium. 

For the case of linear gain ("f a < 0) the wave equation predicts an unlimited 
growth of the field amplitude and, hence, of the energy density. This means that a 
stationary lasing state is not possible in this case and, therefore, the limit S7 — must 
strictly not be taken in Eq. ^ . Such a behavior of linear gain is expected only during 
the exponential intensity growth shortly after the onset of lasing. A complete theory of 
random lasing must, therefore, take into account either the full temporal dynamics of 
the system, or in a stationary state additional surface loss effects must compesate for 
the gain in the medium (see section 3). Nevertheless, we can extract a characteristic 
size of a stationary lasing spot from this theory by requiring that the stationary lasing 
state has been reached locally, i.e. within a finite subvolume of the system: Causality 
requires that the pole of P^(Q, J7), Eq. ([2]), as a function of f2 resides in the lower 
complex £1 half-plane. For -f a < this is possible only if all the diffusive modes allowed 
inside a given lasing spot have a wavenumber Q > Q m in — \J Re(—j a /D). This, in 
turn, requires that the spot size is 



Qmin y/~Re(— "fa/ D) 

It is the characteristic, maximal size of a spatial region over which diffusive modes 
can be causally correlated in the stationary lasing state. We conjecture that, hence, 
this size is to be identified with the lasing spot size observed experimentally [15] [13] 
in diffusive random lasers. Since according to the microscopic theory [30] the growth 
rate (—7a) is, for small linear gain, proportional to the avergage gain in the medium, 
"f a cx Ime(r), we predict the spot size to be inversely proportional to the gain. 

More generally, despite the fact that the linear gain assumption is not suited 
to describe stationary lasing, it can be used to estimate the laser threshold, i.e. the 
critical pump rate for lasing. Amazingly, this is a rather general remark. For example, 
in a simpler system of a single microsphere with gain, it has been shown |32j . that 
the scattering coefficients calculated within linear response lose their causality just at 
the point where the sphere crosses its lasing threshold. Applied to our random laser 
system, this means that the threshold for lasing within a spot of size R s is reached 
when the transport coefficient —7, determined by the pump rate via the microscopic 
transport theory Refs. [201 E], reaches the value given by Eq. © 

In Fig. [T]we show the numerical evaluation of the spot size R s as a function of 
increasing Ime s for typical parameters, as given in the figure caption. The imaginary 
part of the dielectric constant is a measure of external pumping, since the gain is 
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Figure 1. The spot size R s , Eq. J3J in units of the scatterer radius ro, as 
obtained by causality considerations (see text) as a function of the imaginary 
part of the dielectric constant of the scatterers. The parameter values used are 
ej, = 1, Ree s = 10, scatterer filling fraction v = 30%, light frequency uj/loq = 2.5. 
The ligth frequency is loq = 2nc/ro where c is the vacuum speed of light. The 
data in the inset are taken from Ref. }13j and refer to the spot size of the modes. 



given by the population inversion of the laser. Therefore, larger pumping yields higher 
inversion and leads to a larger Ime s . The calculated spot size displays a qualitative 
agreement with the experimental data |13j (spot size vs. pump intensity / threshold 
intensity) shown in the inset. 

3. Tranport Theory of Random Lasing 

As remarked in section II, a stationary lasing state in a homogenaously pumped system 
is possible only if the system is finite, so that surface loss effects can compensate 
the gain in the medium. To avoid the causality problem, we consider here a three- 
dimensional random laser model with a homogeneously pumped, active medium which 
extends infinitely in the (x,y) plane, but has a finite, constant thickness d in the 
z direction. The laser-active material is described by the semi-classical laser rate 
equations, and the light intensity transport by a diffusion equation. In particular, the 
rate equations for a four-level laser are 

(4) 

= ^ ( _ + _LW_M^2„,,„ (5) 
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N 3 , (8) 

where Ni = Ni(r,t), i — 0, 1, 2, 3 are the population number densities of the 
corresponding electron level (i € {I... 4}), N tot is the total number of electrons 
participating in the lasing process, 7^ = are the transition rates from level 

i to j, and 7„ r is the non-radiative decay rate of the laser level 2. 7p = 1/rp is 
the transition rate due to homogeneous, constant, external pumping. Furthermore 
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n p h = Nph/Ntot is the photon number density, normalized to N tot . In the stationary 
limit (i.e. d t Ni = 0), the above system of equations can be solved for the population 
inversion n 2 — N 2 /N tot to yield (732 and 710 assumed to be large compared to all 
other rates) 

n 2 = — — 7 7T, (9) 

IP + Inr + 721 {n p h + 1) 

The photon number density (light intensity), normalized to N totl n p h — N p h/N tot , 
obeys the diffusion equation |23j . 

d t n ph = D V 2 n ph + ^ 2 i{n ph + l)n 2 , (10) 

where the last term on the r.h.s. describes the intensity increase due to stimulated and 
spontaneous emission, as described by the semi-classical laser rate equations. Since 
in the slab geometry ensemble-averaged quantities are translationally invariant in the 
(a;, y) plane, but not along the z direction, a Fourier representation in the (x, y) plane 
in terms of n p h {Q\\ , z), ri2(Q\\,z) is convenient, 

d t n ph = - DoQ^riph + D d 2 n ph (11) 

[d 2 Q\, ^ 
+ 721 J j2^n ph (Q\\-Q\\,z)n 2 {Q\\,z) + ~f2in2 

We now seek the photon density response function P(Q\\, z, fi), which describes the 
response of the photon density, n p h, to the distribution of the population inversion, 
ri2, in order to determine the transport coefficients. In the stationary case {dtn p h = 0) 
and in the long-wavelength limit along the (x,y) plane (Qii — > 0), the z derivative in 
Eq. ([TT]) can be expressed without derivatives in terms of n p h and n 2 only. Plugging 
this back into Eq. (| 1 1 1) yields, 

d t + D Q 2 ■ 721712 



n P h(Q\\\. z ,t) = j2in 2 (Q\\\, z ,t) (12) 

and, hence, after Fourier transform w.r.t. time, the diffusion form of the density 
response function, 

P E (Q ]b z,ty= ' 721 , (13) 

il + tQ^D + i£ Z D 

where from Eq. (fT2l) the correlation length ^ is defined as the real, positive quantity, 



£=J-^- (14) 
y 721 n 2 

As seen from Eq. (fTB"]) the pole structure of Pe in this finite-size, diffusive system is 
perfectly causal. The square of the correlation length £ remains positive, indicating 
an effective loss out of a given Q\\ mode. This is due to the loss of intensity at the 
surfaces. Additionally, the mass term becomes less and less significant as the laser 
intensity in the sample builds up, because the relative population inversion clearly 
obeys n 2 < 1 whereas the relative photon number is not restricted. 

Since for homogeneous pumping the averaged photon density does not depend on 
x or y, Eq. (jll[) simplifies in the stationary limit to 

D d 2 n ph = - J2i{n p h + l)n 2 (15) 
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Figure 2. The following quantities are shown as a function of z for different 
values of the pump rate •yp: a) the photon number, which increases monotonically 
with increasing pump rate and has its maximum in the film center (z = 0); b) 
the population inversion, which is inverse proportional to n p ^, see Eq. 10; c) 
the correlation length (spot size), which clearly behaves non-monotonically with 
increasing pumping. For all panels the diffusion constant is Do = lrf 2 72l- 



and n p h(z) is finally determined via Eq. ([§]) by the regular differential equation, 

~2 / \ 721 (TP/721) na , 

d z n ph {z) = -— (7p/72i) . (16) 

n ph (z) + l 

Eqs. (fT6)l ([9]) and (|T4|) comprise the complete description of the spatial photon density 
profile perpendicular to the lasing film and the intensity correlation length (spot size) 
parallel to the film. 

Numerical evaluations of Eqs. ([To]) . (TH| and © are shown in Figs. and 
[3] In Fig. [2] the photon number n p h(z), population inversion n 2 (z) and correlation 
length £(z) are shown as a function of z for different values of external pumping, 
characterized by the pumping rate jp. The value of the diffusion constant was chosen 
to be Dq — Id 2 721, where d is the width of the film. In panel a) of Fig. [2] the photon 
number displays a monotonically increasing behavior with increasing pumping. The 
maximum of the intensity resides in the center of film {z = 0), since this is the 
position farthest from the boundaries, and therefore with lowest loss of intensity. The 
population inversion, Eq. ([9]), behaves inverse to n p h(z) 1 see Eq. (j9]). In contrast to 
this rather expected behavior, the correlation length £(z) as given by Eq. (fT4|) exhibits 
a non-monotonic behavior with increasing pumping. For pumping rates "fp < 721 the 
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Figure 3. The figure shows the following quantities as a function of the pumping 
rate P at the film surface (z = ±0.5d): a) the photon number which displaying 
saturation for strong pumping; b) the population inversion, also saturating; c) 
the correlation length (spot size) showing non-monotonic behavior. The inset of 
panel c) displays the 1 / y/*yp dependence of the intensity correlation length £ on 
the pump rate above threshold (-yp > 1) (sec discussion in the text). 



correlation length increases but for pumping rates 7p > 721, £ is decreasing. The 
equality between jp and 721 marks the situation where electrons are as fast excited 
into the upper laser level as they relaxate to lower levels. Therefore this characterizes 
the lasing threshold. Available experimental data [HI [13] also report a decreasing 
behavior of the spot size above threshold. Measurements of the intensity correlation 
length below threshold have not yet been reported. 

The same quantities are shown Fig. H as a function of external pumping 7p 
at the surface of the random laser. Photon number and population inversion both 
display saturation behavior. Panel c) of Fig. however, exhibits the non-monotonic 
behavior of the correlation length. This plot is to be directly compared to experimental 
data P2], as e.g. shown in the inset of Fig. [5] There is a good qualitative and even 
quantitative agreement between calculated and measured spot size. The inset of panel 
c) exhibits shows that the dependence of the spot size £ on the pump rate 7p above 
threshold (7p) is predicted to be 

£( 7 p) = £(oo) + a/V7F (17) 
with a proportionality constant a. This result is open for further experimental tests. 
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4. Conclusion 

We have discussed how a linear response theory for light transport, including self- 
interference effects, in disordered media with linear gain predicts threshold behavior 
of the intensity. Even more interestingly it also predicts a characteritic, average radius 
of lasing modes, dictated by causality. We identify this length scale with the spot size 
of the random laser as measured in experiments [15l H3] and find qualitatively good 
agreement. Further, we have proposed an analytical transport theory for random 
lasing in finite systems. A finite system size is necessary for surface loss to compensate 
the gain in the medium and, hence, to stabilize a stationary lasing state. In particular, 
we consider a slab geometry where in the medium the light intensity propagates 
diffusively and the loss through the surfaces is included by appropriate boundary 
conditions. The theory allows for the first time for an analytical calculation of the 
intensity correlation length of this system, describing the spatial extent of a mode 
(spot size). The spot size is predicted to behave non-monotonously as a function of 
external pumping, i.e. increasing below and decreasing above the laser threshold. A 
comparison with experiments reveals qualitativly good agreement. Our prediction of 
its functional dependence on the pump rate is open to further experimental tests. 

Our future work will include solving the semi-analytical light transport theory 
with self-interference contributions when the system is self-consistently coupled to the 
laser rate equations. 

Useful discussions with H. Cao, P. Henseler, B. Shapiro, and C. M. Soukoulis 
are gratefully acknowledged. This work was supported in part by the Deutsche 
Forschungsgemeinschaft through grant no. KR1726/3 (R.F., J.K.) and FG 557. 
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